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Abstract 

We argue that there is no essential violation of universality in the 
continuum limit of mixed and 0(n) lattice sigma models in 2 

dimensions, contrary to opposite claims in the literature. 
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1 Introduction 


In this paper we consider two-dimensional mixed isovector-isotensor 0(n) 
sigma models described by a lattice action of the kind 

^(5) =PyY,{l- s.s,+^) + ]-f5T E (l - (SxS,+^)2) , (1) 

with = 1. The sums run over the nearest neighbor sites. This provides 
a possible lattice discretization for the continuum 0(n) non-linear sigma 
model, 

^cont ^ y ^2^ (a^S(x))2 (2) 

with /3 = /3v + Pt- 

According to conventional wisdom, different lattice regularizations (pre¬ 
serving the crucial symmetries) yield the same continuum field theory (“uni¬ 
versality”). For the case of the action (|I|), Caracciolo, Edwards, Pelissetto 
and Sokal ii, however, question this assumption and in particular state 
that the pure sigma model (/3 t = 0) and the pure model (/3y = 0) 

have different continuum limits for (5 ^ oo. Since the notion of universality 
plays an essential role in the theory of critical phenomena it is worthwhile 
to consider this question again. In this paper we will explain how the pecu¬ 
liar features observed in the model (|l|) can be understood in the framework 
of the conventional picture. We wish to stress, however, that our scenario 
is (for the most part) based on plausibility arguments, for which rigorous 
proofs are unfortunately still lacking. 

A related problem concerns the mixed fundamental-adjoint action in 
pure SU(n) gauge theory in 4 dimensions. The generally accepted belief 
is that there is a universal continuum limit for these theories. However, we 
shall not discuss this model here. 

The paper is organized as follows. In section 2 we consider a class of pure 
Rpn-i \Yg £j.g^ describe some general properties and then go on to 

discuss the continuum limit. Section 3 presents an investigation of perturbed 
Rpn-i 2 ^Q(^g}g^ paying special attention to their expected continuum limit. 
In particular, we argue there is no contradiction to the general understanding 
of universality. Finally in section 4 we outline some calculations supporting 
our general scenario. 
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2 The RP^ ^ models 


2.1 Some general properties 

The standard action of the model is 

^t(5) = ^/3E(i-(S-W)- (3) 

It has, compared with the 0(n) model, an extra local Z 2 symmetry: it is 
invariant under the transformation 

Sx gxSx, where gx = ±1. (4) 

As a consequence, only those quantities have non-zero expectation values 
which are invariant under this local transformation. In particular the isovec¬ 
tor correlation function vanishes; 

(SxSy) = 0 for X ^y. (5) 

The simplest local operator with non-vanishing correlation function is the 
tensor 

{TfTf) + 0. (6) 


This behavior seems completely different from that of the 0(n) sigma 
model, so that one might expect drastic differences in the physics described 
by the models. This is indeed true for the theories with finite lattice spacing, 
but below we shall argue that in the continuum limit this difference becomes 
insignificant, and can be resolved by consideration of nonlocal variables. 


2.2 Defects and phase structure 

For convenience, we introduce the notation Uxy = SxSy for the scalar prod¬ 
uct of two spins. Further, for any path V on the lattice define the observable 

WiV)= n (7) 

<x,y>G'P 
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where < x,y > denotes the link joining two neighboring points x and y. 

Consider a configuration of the model. One says that it has a 

defect associated with a plaquette p (or a site on the dual lattice) if 

W{dp) < 0 , (8) 

where dp is the boundary of the plaquette. The defects are endpoints of 
paths on the dual lattice formed by those dual links with Uxy < 0, where 
X, y are the two sites on the corresponding link. Due to the local gauge 
invariance, only the position of the defects is physical, while the paths can 
be moved by a gauge transformation. 

Like the vortices in the two-dimensional XY model Q , these defects play 
an essential role in determining the phase structure of the model at 

finite (3 Q. Some of these aspects are discussed by Kunz and Zumbach |^. 
The activation energy of a pair of defects grows logarithmically with their 
separation r. The standard energy-entropy argument Q then predicts a 
phase transition at some finite Pc- For P < Pc the defects are deconfined, 
while for P > Pc they appear in closely bound pairs. This difference is 
expected to show up in an area or perimeter law (for P < Pc and P > Pc 
respectively) of the “Wilson loop” expectation value (fT(£)) for large loops 

^ i- 

We see this in a large n limit of the RP""”^ model |^, |8[. There the phase 
transition is demonstrated to be first order. Furthermore, one verifies the 
expected “Wilson loop” signal: in the leading order, (VF(£)) = 0 for /3 < Pc, 
while {W{C)) = exp{—7(/3)|£|} for P > Pc, with |£| the perimeter of C. 

For finite n, however, the situation is not at all clear. The discussion of 
the nature of the critical point at finite P has a long history I, 0,0,ii. 
All MC simulations show that approaching Pc from below the correlation 
length starts to grow drastically. However, the various authors disagree 
concerning the nature of the transition, the variety of opinions based merely 
on theoretical expectations (and prejudices). We shall return to this question 
later. 

In the following we will discuss the possible continuum limits. We shall 
argue that at finite P the correlation length in the RP”'“^ model always 
stays finite, and the critical point at /3 = oo is equivalent to that of the 0(n) 
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model. 


2.3 Equivalence of the RP” ^ and 0(n) models in the contin¬ 
uum limit 


Consider a more general form of the lattice RP"' ^ action: 

AT{S) = f3 Y. f{u,y), (9) 

<x,y> 

where the function f{u) satisfies the following properties: 

f{-u) = fiu), /( 1 ) = 0 , /'( 1 ) = -1 ( 10 ) 

and f{u) is monotonically decreasing for 0 < u < 1. We assume a weaker 
form of universality: any of these choices yields the same continuum limit 
as /? —> oo. (Actually, even less will be sufficient — one can keep hxed the 
form of f{u) for uq < |m| < 1 to be the standard one.) 


Let us now introduce a chemical potential /r of the defects modifying the 
Boltzmann factor by exp(—/rn^e/) where Udef is the number of defects. At 
/i > 0 the defects are suppressed and at /r = oo no defects are allowed. 


Take first the fi = oo case. As was done by Patrascioiu and Seiler |R], 
one can dehne Ising variables e^; = ±1 by 


= sign{VP {Vxox)} , (H) 

starting from a fixed site xq and going along some path Vxqx connecting xq 
to X. Due to the absence of defects, ex does not depend on the path chosen. 
For two nearest neighbor sites one has 


^X^X+fJ, Sign('U2;3;+^) . (l^) 

Introduce now a new 0(n) vector 

(T X — (l^) 

This has the property that UxCrx+y = > 0 for nearest neighbors. 

The dynamics of the (Tx held is described by the modified 0(n) action 

Av{cr) = i^xCTx+y,) (14) 

x,ii 
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with 


(15) 


= for«>0, 

I oo for u < 0 . 

We also assume that the continuum limit {(3 ^ oo) for this action is the 
same as for the standard 0(n) action (universality within the 0(n) model). 

The model described by (|9|) at /i = oo and the corresponding 0(n) 

model given by ( p^ ) are equivalent in the continuum limit in the following 
sense: all gauge invariant quantities (such as the tensor correlation function 
or a Wilson loop of scalar products) in the model are exactly the 

same as in the 0(n) model, while all non-gauge invariant quantities vanish 
in the RP”“^ model. In particular, for the vector correlation function 

(SxSy) = {exey){axcry) = 0 ioi X ^ y (16) 

since {ex^y) = hxy The vector of the RP"""^ model can be thought of 
as a product of two independent fields, the “true vector” ax and the Ising 
variable e^,; one is described by the corresponding 0(n) model, while the 
other by an Ising model at infinite temperature. 

We return now to the case of RP""”^ model at finite y. With increasing y 
the average defect density is decreased. Defects tend to disorder the system, 
therefore it is very plausible to assume that the correlation length (in the 
tensor channel) grows with increasing y. Since at // = oo the RP**”^ model 
is equivalent to the corresponding 0(re) model at the same (5, one concludes 
that the correlation length at /x = 0 is bounded by that of the 0(n) model. 

Assuming further that, according to the standard scenario, the 0(n) 
model has a finite correlation length for finite (3, it follows that the RP”'”^ 
model cannot have a phase transition (at finite f3) with diverging correlation 
length. 

The latter is in agreement with the large n result |8| mentioned above, 
which predicts a first order transition. The explanation for the seemingly 
divergent correlation length observed in MC simulations could be the follow¬ 
ing. For (3 < (3c the defects strongly disorder the system and cause a small 
correlation length. Above (3c, however, the role of the defects decreases 
rapidly with increasing p. As the defects become unimportant the correla¬ 
tion length approaches that of the 0(n) model. The numerical simulation 
of the RP^ model gave Pc = 5.58 which in the 0(3) model corresponds to 
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a correlation length ^ ~ 10^®! A sharp transition or a jump to a huge value 
is therefore is not unexpected. This transition is, however, associated with 
the non-universal dynamics of the defects, not with the universal continuum 
limit of the theory. 

To establish the equivalence of the model (at fi = 0) with the 

0(n) model in the continuum limit it suffices to show that the defects do 
not play any role in the (3 ^ oo limit. The defects (or rather pairs of defects) 
have finite activation energy which depends on the distance r between the 
two defects as const + ^vrlnr. The constant contribution coming from the 
neighborhood of the defects depends strongly on the actual form of the 
function f{u) in (P), more precisely on the values of f{u) for small |u|, say| 

< 0.5. Because the defect pairs have finite activation energy Eq, they 
are exponentially suppressed by exp(—/3illo). The subtlety here is that the 
correlation volume, oc exp(47r/3) (for n = 3), is also exponentially 

large, and pairs of defects with limited relative distances will occur in this 
volume if their Eq is small enough^]. These could be, however, considered 
as local — i.e. non-topological excitations on the scale of ^(/9), and we do 
not expect that they significantly influence the P ^ oo limit. The argument 
becomes even simpler if one changes the form of the action by pushing up 
the values of f{u) for < 0.5 to have Eq > dvr for all defects. In this case 
the defects are practically absent in the whole correlation volume |^. 

As a concrete realization of the modified model we take 

fiu) = ^(1 -u^) + q- max(uo - 0). (17) 

Here q > 0 and we choose Uq = 0.8 for definiteness. A simple numerical 
investigation shows that for g = 10 the activation energy for neighboring 
defects is Eq « dvr. (Of course, nothing forbids taking q = oo — it will still 
define the same continuum theory.) 

By similar modifications of the action it might well be possible to bring 
the correlation length down to reasonable values, so that the phase diagram 
could be reliably investigated numerically (also in the mixed RP”'“^/0(n) 

® It is easy to show that around a defect at least one of the four links has < 0.5. 

^For the standard RP"'“^ action the minimal activation energy is = 2.14 

® Obviously this argument does not apply if the correlation length becomes infinite 
already at Hnite /3 as suggested in ref. |12|]. 
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model). This would imply that the huge correlation length around the point 
where the defects start to condensate for the standard model is rather 

“accidental”. 


3 The perturbed RP” ^ model 


Consider the perturbed RP” ^ model 

A{S) = (5TY.f (SxS,+^) + I3vY.9 (SxS,+^) (18) 

in the limit (3t oo, j3v = fixed. Here f{u) satisfies (0), while the pertur¬ 
bation g{u) can, without loss of generality, be taken to be odd: 

g{-u) = -g{u). (19) 

The action (|l]) is, of course, (up to an irrelevant constant) a special case. At 
Pt ^ oo the scalar product Sa;Sa;+^ is forced to be around -|-1 or —1, i.e. 
1 - ( S , S ,+^)2 = 0 { 1 / Pt ). 


Let us now assume that Pt is large enough or the form of /(u) is chosen 
such that the defects are completely negligible (as in the example of ( 0 ) 
for q > 10). For configurations with no defects one can introduce the Ising 
variables in a unique way and define the “true vector” field ax as in ([T^. 
Separating the sign of g{u) by 


g{u) = -sign{u)goi\u\) = -sign(n) [50(1) + 5 'o(l)(kl - 1 ) + ■ • •] , 
we obtain 

>1(S) = Av{(7) + Aisingie) + Aint{e, a), 

where 

Av(o-) = Pj^fv(o-xO-x+g), 

X,/l 

A/singi^) — J ^ ^ 

X,fM 

Ainti^j ^') — ^ ' (-x^-x+g, [ 5o(^)(^ T . . .] . 

x,g 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 
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Here (3 = Pt, J = Pv9q{^) and fv{u) is as in (|I5|). Note 1 — axC^x+ij. = 
0{\/P) and hence the interaction term Aint{^,o') goes effectively to zero as 
P oo. 

Consider first the simple case when g{u) = —sign(u), i.e. go{u) = 1. In 
this case the two systems decouple exactly while the specific behavior of the 
vector and tensor correlation functions still persists. Since the correlator 
(SxSy) factorizes: 

(S^Sy) — {Ax^y) 1 

for J < Jc one has 

ms = me + mo- and mr = 2mo-, (26) 

where the masses are defined through the exponential decay of the corre¬ 
sponding correlators. Although the tensor mass is smaller than twice the 
vector mass, m^ < 2ms, one can not conclude from this that there is a pole 
in the tensor channel (in contrast to the pure 0(n) model), as suggested in 
ref. 1^. Since both m^iP) and mPJ) go to zero as P and J approach their 
critical values, the ratio 


mx 2ma 

fns + ITT'e 


(27) 


can be fixed at any value r G [0, 2] by properly approaching the point (Jc, oo) 
in the (J, P) plane. 


For J > Jc the Ising field ex develops a non-zero expectation value hence 
in this case ms = m^ and mx/ms = 2. Note that for finite p the phase 
transition around J = Jc is observed only in the non-local variable tx not in 
the original variable S^; whose correlation length remains finite ai J = Jc- 

Following the argument in refs. ii one would conclude that around 
the point (J,/?) = (Jc,oo) one could define seemingly inequivalent theories 
differing in the ratio mx/ms 0. Although this is formally true, the corre¬ 
sponding theory is neither really new nor interesting. In particular, all the 
tensor correlation functions are the same as those in the corresponding pure 
0(n) model. 

®The masses measured in Q are not the true masses, but those defined through the 
second moments; it is however generally believed that the qualitative picture remains 
unaltered. 
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With the choice g{u) = —n, i.e. go{u) = |tt| (as in [^) the situation 
is more complicated since there is an interaction between the two systems. 
However, as mentioned above, the effective strength of the interaction goes 
to zero as /? —> oo, hence it might well happen that in the continuum limit 
one recovers the previous situation. 

Note that the presence or absence of the interaction is not connected 
with the behavior of g{u) around u = +1 (which is responsible for the 
0(n) continuum limit f3v oo) but rather with the difference in behavior 
around u = +1 and u = —1. For example, g{u) = |(1 — n^) + c9{—u) (not 
antisymmetrized in this case) where c > 0 and 6 is the step function, is a 
perfectly acceptable discretization of the 0(n) model for /3y ^ oo and it 
produces no interaction, Aint = 0. On the other hand, g(u) could be chosen 
to have, say, a local maximum at u = +1 instead of a minimum, which 
would completely destroy the /9y —> oo behavior but would still have the 
same interaction pattern as for the case g(u) = —u. 

In this sense, the phenomenon around the point (Jc, oo) is the conse¬ 
quence of perturbing the model by a term breaking the local Z 2 

symmetry, rather than its mixing with the 0(n) model. 


4 Some analytic studies of the mixed model 

Let us set /3y = (1 — uj)n/f and (St = ojnjf for the bare couplings in (Q). 
There are various analytic studies which shed some light on the physics of 
this model. Among these are the ordinary perturbation theory / —> 0 and 
the 1/n approximation. 


4.1 Bare perturbation theory 

One interesting exercise is to compute the spectrum for a finite spatial extent 
L. For the tensor mass rriT to second order in bare perturbation theory, one 
finds 

mT{L)L = / + /2i{(n - 2)R{L/a) + [1 + u^{n + l)]P{L/a)] + 0{f) (28) 
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and to this order the vector mass my is given by 


mv{L) = mT{L). 


( 29 ) 


In (25) the functions R,P are given by finite sums over lattice momenta. 


The relation (29) holds before the continuum limit has been taken (there 
are no lattice artifacts in the ratio to this order). Furthermore, the ratio is 
independent of cn, which is certainly consistent with notions of universality 
(the continuum limit is taken here in finite volumes). The ratio (|2^ has been 
shown to hold in the 0(n) model for small volumes, in the continuum limit 
to third order in the renormalized coupling by Floratos and Petcher fl^ . 
Indeed there, to this order, the mass of the tensor of rank k is proportional 
to the eigenvalue of the square Casimir operator: 


ruk = Mk{n + k — 2) 


(30) 


with M independent of k. In finite volumes the spectrum is discrete and 
there is a finite gap between rriT and 2my; this gap is expected to close 
as L —> oo where a cut develops starting at 2my. We have numerically 
computed the mass of the tensor as well as that of the “true vector” in the 
Rpn-i g^g defined in sect. 2, in small volumes; the results agreed well 

with the above formulae. 


One can also use (pSD to determine the ratio of A-parameters. For this, 
it suffices to know the continuum limit (a/L —> 0) behavior of R, P: 


R{L/a) ~ 


^|ln(L/a) - ln(7r/\/2) + 7 e|, 


(31) 


P{L/a) ~ - 


(32) 


with 'jE Euler’s constant. Denoting A(a;) the lattice A-parameter for a model 
with given w, 


A(uj) ( ujTrin + 1) 'i 


(33) 


follows, in agreement with the result in ref. H- 


Comparing the two theories in infinite volume, Caracciolo and Pelissetto 
[jl^ also found that the RP””^ and the 0(n) models have (apart from the 
redefinition of the coupling) the same perturbative expansion. 
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4.2 1/n Expansion 


The 1/n expansion for the mixed model was to our knowledge first inves¬ 
tigated by Magnoli and Ravanini ||^]. We disagree, however, with some of 
their final conclusions. To discuss this, we first introduce a few formulae. 
After introducing auxiliary fields A^(x),t(x) to make the integral quadratic 
in the spin-fields and then performing the Gaussian integral, the partition 
function in the absence of external fields, takes the form 

Z = const- f J|dA^(x) J|dt(x)exp|-^5eff|, (34) 


with the effective action 


5'eff 



-|- it{x)] -i- tr In A4, 


(35) 


where M is the operator 

M = s + it + ^ I-9* -I- u![And*^dfj, - (cl*A^)(l - 9*) -|- A^]|. (36) 

Here 5/^(9^) denote the lattice forward (backward) derivatives. One first 
seeks a stationary point of at constant field configurations A^{x) = 1 — 6, 
t{x) = const. Demanding a saddle point at t = 0 gives a relation for the 
constant s in ( |35[ ) as a function of 6. With s fixed in this way, one seeks 
minima of S'efr as a function of 6. 

For a; = 1 (the pure model), the extremal points are shown in fig. 

1. In this case there is a symmetry 6 —> — 6 [|. Further 6 = 0 is an extremal 
point for all /. For / < 1, the points 6 = 0 are maxima and the non-zero 
values are minima. For / = I"*", 6 = 0 becomes a local (but not absolute) 
minimum and two new local maxima develop. At / = /c(l) ~ 1.046 the 
three minima become degenerate, and for / > /c(l) the minimum at 6 = 0 is 
the absolute minimum. One finds (in the leading order of the 1/n expansion) 
that at this point the tensor correlation length does not go to infinity: there 
is a jump in the order parameter and the phase transition is thus first order. 

^Actually, for the pure case there are degenerate minima, due to the 

local Z 2 symmetry. Elitzur’s theorem is not violated by this approximation — the local 
symmetry is not broken spontaneously. 
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Figure 1: The order parameter b in the 1/n expansion as a function of the 
coupling / for uj = 1. There is a jump at /c(l) = 1.046 from a finite value 
to 6 = 0 shown by the dotted line. 

For a; < 1, the 6 —> — 6 symmetry is broken and the local minimum 
with 6 > 0 is the lowest. For uj only slightly less than 1, the situation is 
as in fig. 2. Here again, at some / = fc{oj) the parameter h undergoes a 
finite jump. There is, however, a critical value of a; = Wc ~ 0.985 below 
which the “S-structure” in fig. 2 dissolves and there is only one extremal 
point for b > 0 for all values of /. In the to — f plane there is thus a 
first order transition line which starts at (1, /c(l)), extends only a little way 
in the plane and ends at a critical point C = {ujc, fd^^c) ~ 1.075). At 
C the vector and tensor correlation lengths remain finite. The transition 
at C is, however, second order since the specihc heat diverges. The cause 
of this in the leading order of the 1/n expansion can be traced back to 
a development of a singularity in the inverse propagator of the auxiliary 
fluctuating t-field|^ at zero momentum at the critical point. The singularity 
in the t— propagator seems to remain for higher orders as well. An infinite 

®Note that the A— and t— fields mix and it is necessary to diagonalize. 
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correlation length in the energy fluctuations does not contradict a finite 
correlation length in the vector and tensor channels; in particular, there is 
no conflict with correlation inequalities. These inequalities state that by 
increasing a ferromagnetic coupling the system becomes more ordered and 
the correlation between any spins increases. Although this assumption looks 
physically quite obvious, it has not been proven rigorously. The increase 
of the correlation function, however, implies the growing of a correlation 
length with increasing ferromagnetic coupling, only when the corresponding 
quantity has a vanishing expectation value. 



Figure 2: The order parameter 6 as a function of / for uj = 0.999. It still 
has a finite jump indicated by dotted line. At a; > cjc = 0.985 the S-shape 
dissolves thus the phase transition disappears. 

Thus, a diverging vector (or tensor) correlation length at the endpoint C 
would contradict a finite correlation length for large (but finite) (3v (asymp¬ 
totic freedom) — on the other hand, a diverging specific heat at C is not 
excluded by these considerations. The above scenario disagrees with that 
of Magnoli and Ravanini Q who argue (based on correlation inequalities) 
that the second order phase transition at the point C is only an artifact of 
the 1/n approximation. 
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Caracciolo, Pelissetto and Sokal [1^| also discuss the (3t/N^ j3v fixed, 
—> oo limit. They obtain a result which is equivalent to eq. above 
(although their interpretation is different from ours). 


In conclusion, it is plausible that the phase diagram, also for finite n, 
is the “standard” one shown in fig. |3|. There is a first order transition 
line starting at the point A of the (5t axis. It ends at the point C where 
the specific heat becomes infinite, but the vector and tensor correlation 
lengths remain finite. In this figure we also indicate the Ising critical point 
B discussed in Section 3. The dotted line starting at point B is the critical 
line of the underlying Ising variable e. This criticality, however, does not 
show up in the correlation functions of the original variable S. 



Figure 3: The phase diagram for the mixed RP" ^ - 0(n) model. 
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